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INTRODUCTION 

Despite  the  unfortunate  circumstances  of  being  unable  to  use 
the  Vapor  Phase  Synthesis  Laboratory  for  the  months  of  October  through 
January  because  of  construction  in  that  portion  of  the  McCullough 
building,  considerable  progress  has  been  achieved  during  the  past  six 
months..  This  report  focuses  on  five  individual  work  areas  (i)  synthesis 
of  SiO^,  (ii)  synthesis  of  SiC,  (iii)  synthesis  of  Pd2Si,  (iv)  electri¬ 
cal  diagnostics  of  films,  and  (v)  computer  simulation  of  film  synthesis. 
In  the  Si02  area  we  have  found  that  the  film  thickness  is  proportional 
to  T-1  with  an  activation  energy  of  1.5  kcal/mole,  that  the  film  uni¬ 
formity  increased  with  P  increases  and  that  IR  transmission  spectra 

°2 

showed  a  spectral  shift  to  higher  wave  numbers  as  the  deposition 
temperature  increased.  In  the  SiC  area,  single  source  deposition  led 
to  6- SiC  films  with  a  (111)  texture  and  measurements  led  to  data  on 
(a)  film  stoichiometry,  (b)  refractive  index,  (c)  film  structure  and 
morphology,  and  (d)  sheet  resistivity.  Steady  progress  has  been  made 
and  the  addition  of  tin  as  a  mobility  enhancer  has  produced  the  best 
quality  film  to  date.  Work  has  begun  on  the  synthesis  of  (Pd,  Pt,  Ir) 
silicides  and  dual  source  deposition  of  Pd  and  Si  has  developed  the 
composition  and  phase  spread  in  that  binary  system.  In  the  area  of 
electrical  diagnostics,  theoretical  review  of  the  I-V,  C-V  and  AC 
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conductivity  characteristics  for  p/n  junction,  hetero junction  and  metal 
insulator  semiconductor  diodes  has  been  conducted  and  experimental 
measurements  has  been  made  on  SiOx  MIS  diodes.  Another  new  area  of 
study  has  been  computer  simulation  of  film  synthesis  where  a  number  of 
new  software  programs  have  been  developed.  In  addition,  a  special 
parametric  study  of  a  multilayer  structure  representing  the  Cu/Ag 
system  has  been  conducted  and  special  potential  energy  functions  have 
been  developed  for  the  Si-0  system  in  order  to  study  its  unique 
structural  features. 

A.  EXPERIMENTAL  RESULTS 

1 .  Reactive  Sputter  Deposition  of  Oxide  Films 


A  number  of  experimental  modifications  have  been  made 
relative  to  the  sputtering  chamber  configuration  reported  on  previ¬ 
ously^  (see  Fig.  1).  An  Argon  ion  gun  was  added  for  the  removal 
of  native  oxide  prior  to  deposition.  This  necessitated  making 
flexible  power  leads  for  the  substrate  heater  so  that  it  could  be 
rotated  from  under  the  ion  gun  to  a  position  beneath  the  target. 

Nickel  heat  shields,  having  relatively  inert  oxides  at  high 
temperatures,  were  fabricated  for  the  heater.  In  addition,  a 
rectangular  quartz  oxygen  dispersion  ring  was  made  in  order  to  enhance 
deposition  uniformity  over  the  rectangular  area  dictated  by  the 
heater  (see  Fig.  1). 


To  ensure  native  oxide  removal,  an  investigation  of  ion  gun 

operating  conditions  was  made.  The  critical  variables  were  determined 

to  be  (a)  ion  gun  to  substrate  distance,  (b)  ion  gun  power,  and  to  a 

lesser  extent,  (c)  chamber  pressure.  The  thermal  oxide  etching  rates 

for  various  ion  gun  etching  parameters  are  given  in  Table  1.  Etch  rates 

of  approximately  14  X/min.  were  obtained  at  a  4  inch  operating  distance 

for  the  ion  gun  at  2kV  and  2mA  and  the  Argon  chamber  pressure  at  4vun 

and  8ym.  At  a  5  inch  operating  distance,  2kV  and  2mA,  the  etch  rate 

was  only  3  X/min.  for  an  Argon  chamber  pressure  of  7.25ym.  At  a  3  inch 

© 

operating  distance,  etch  rates  up  to  22  A/min.  were  obtained. 

The  thermal  oxide  etch  profile  was  examined  ellipsometrically 
and  is  shown  in  Fig.  2.  Significant  etching  was  accomplished  over  a 
region  ^1  inch  diameter.  For  long  etches,  the  profile  was  quite 
steeply  sloped  with  a  spot  'v  one- fourth  inch  diameter  exhibiting  a 
maximum  and  uniform  etch.  The  effect  of  Argon  chamber  pressure  on 
etch  spot  size  or  etch  rate  was  negligible  for  Argon  pressure  variations 
in  the  range  4ym-8ym.  It  was  also  noticed  that  the  film  index  of 
refraction  (n^)  varied  with  the  amount  of  material  etched  (see  Fig.  3) 
which  indicates  possible  damage  to  the  thermal  oxide  (and  also  the 
ellipsometric  sensitivity  to  such  changes). 

A  series  of  oxide  depositions  were  made  in  order  to  investigate 
the  influence  of  substrate  temperature  as  a  deposition  parameter  with 
the  other  conditions  such  that  SiO^  should  form.  These  deposition 
conditions  are  listed  in  Table  2.  The  temperature  at  the  substrate 
was  varied  from  25°C  to  660°C  using  a  single  element  quartz  heater 


and  the  resulting  oxide  thicknesses  showed  a  non-linear  decrease  with 
increasing  temperature .  X-ray  diffraction  indicated  that  all  the  films 
were  amorphous  and  electron  microprobe  analysis  showed  the  oxygen/ 
silicon  ratio  to  be  about  2/1  except  in  the  case  of  run  #81-218  where 
the  index  of  refraction  was  1.36  (indicating  porosity  or  excess  oxygen) 
and  the  0/Si  ratio  was  ^5/4.  Since  the  microprobe  samples  a  depth  of 
'''O.Sy,  it  is  likely  that  the  Si  substrate  was  involved  in  the  results 
of  samples  #81-220  and  #81-221. 

Because  nickel  has  a  high  sputtering  yield,  the  silicon  mask 
shown  in  Fig.  4  was  fabricated.  A  piece  of  quartz  used  to  define  a 
thickness  step  for  step  profile  thickness  measurements  (ellipsometric 
calibrated) .  Using  the  ion  gun  parameters  reported  above,  conditions 
of  operation  to  achieve  etching  were  chosen  (3  inch  gun  to  substrate 
distance,  2kV/2mA  beam  power,  4.25ym  Argon  pressure,  10  minute  etch) 
and  kept  constant  for  the  next  set  of  depositions.  A  method  for 
accurately  locating  the  ion  beam  center  was  developed  which,  coupled 
with  a  knowledge  of  the  etch  profile,  allowed  localized  areas  of  each 
substrate  to  be  compared;  i.e.,  native  oxide  removed  areas  with  native 
oxide  present  areas. 

Following  this,  a  second  series  of  depositions  were  made  at 
varying  substrate  temperatures  concentrating  on  the  region  of  apparent 
non-linear  behavior  mentioned  above  (100°C-200°C) .  In  this  case, 
comparing  regions  of  maximum  oxide  thickness  on  each  substrate,  the 
linear  dependence  of  film  thickness  vs.  T~*  shown  in  Fig.  5  was  observed 
The  activation  energy  for  this  process  is  ^1.5  kcal/mole.  It  is 
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believed  that  film  #81-218  (240°C)  is  anomalously  too  thick,  perhaps  due 
to  an  error  in  experimental  procedure.  Again  X-ray  diffraction  indicated 
that  all  the  films  were  amorphous. 

To  investigate  whether  or  not  the  oxide  uniformity  could  be 

improved,  the  oxygen  flow  rate  was  increased  leading  to  the  result 

indicated  in  Fig.  6.  At  the  higher  oxygen  flow  rate,  an  improvement 

in  the  index  of  refraction  constancy  over  the  substrate  surface  was 

noted;  i.e.,  from  r|£  =  1.45  ±  0.07  at  PQ  =  0.38pm  to  *  1.458  ±0.003 

at  P  *  0.53pm.  Unfortunately,  the  Si  target  was  oxidized  at  the  higher 
°2 

PQ  resulting  in  a  50%  drop  in  oxide  thickness  for  equal  deposition 
times.  Just  recently,  a  new  Argon  outlet  became  available  which  will 
permit  the  release  of  Argon  radially  around  the  surface  of  the  Si 
target.  This  should  provide  much  better  isolation  of  the  target  from 
oxygen  released  at  the  substrate  and  hence  permit  higher  PQ  values  to 
be  used  without  target  oxidation. 

Auger  studies  have  shown  the  earlier  SiOx  films  to  be  uniformly 
dispersed  on  an  atomic  scale  rather  than  dispersions  of  Si  clusters  in 
Sit^  such  as  found  in  SIPOS  material.  Infrared  transmission  spectra 
studies  on  the  amorphous  SiC^  films  show  that  those  deposited  at  higher 
temperature  (T  'WOO'C)  have  a  shifted  spectra  (to  higher  wave  numbers) 
compared  to  those  deposited  at  lower  temperatures  (T'v25°C).  A  typical 
IR  spectrum  for  amorphous  SiO^  is  shown  in  Fig.  7  while  the  data  compari¬ 
son  for  various  T  is  shown  in  Fig.  8.  Thus,  we  expect  that  the  relative 
perfection  of  these  films  can  be  simply  monitored  by  following  the  IR 
transmission  peak  shifts. 


In  a  previous  report, ^  Table  III  describes  the  synthesis 
conditions  for  SiOx  compounds  by  atomic  layer  growth.  Here,  the  sub¬ 
strates  are  mounted  on  a  turntable  rotating  at  a  fixed  rate  controlled 
to  0.1%  and  first  passed  beneath  a  sputter  source  of  the  metal  atom 
before  reaching  a  region  where  the  newly  deposited  layer  was  exposed 
to  the  reactive  gas,  etc.  The  individual  metal  layer  thickness  was 
controlled  by  the  rotation  rate  of  the  substrate  table  and  the  sputtering 
power  of  the  metal  target.  The  partial  pressure  of  the  reactive  gas 
over  the  substrate  coupled  with  the  rate  of  substrate  motion  determined 
the  exposure  of  the  deposition  surface.  For  this  reason,  we  can  con¬ 
sider  this  a  cyclic  deposition/reaction  mode  of  film  formation. 

The  process  data  for  three  of  these  samples  are  listed  in 
Table  3  while  Figs.  9  and  10  present  measurements  on  film  thickness, 
tp  and  index  of  refraction,  r|£»  as  a  function  of  radial  distance  on 
the  wafer.  The  film  thickness  was  determined  by  the  Tally-step  method. 

As  expected,  the  film  thickness  increases  and  the  index  of  refraction 
decreases  as  the  rotation  period  decreases  (more  complete  oxidation) . 

The  remarkably  uniform  index  of  refraction  observed  at  the  fastest 
rotation  rate  is  of  significant  interest. 
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In  Table  IV  of  a  previous  report, v  J  the  synthesis  process 
parameters  for  several  ZnOx  films  were  tabulated  and  the  electrical 
resistance  of  one  of  them  (#024  )  was  presented  (Fig.  6) .  J 


Electron  microprobe  results  on  this  film  as  a  function  of  radial 
position  from  the  center  of  the  wafer  is  given  in  Fig.  11  so  we  see 
that  the  film  composition  changes  by  W5%  from  the  center  to  the 
periphery  while  the  electrical  resistance  changes  by  a  factor  W  over 
the  same  distance.  It  is  the  index  of  refraction  data  of  Fig.  12  for 
this  film  that  is  most  interesting.  It  reveals  the  presence  of  several 
singularities  along  the  radial  coordinate  indicating  that  the  film  is 
rich  in  structural  variations  that  need  to  be  carefully  explored. 

2.  Single  Source  Sputter  Deposition  of  SiC  Films 

In  a  previous  report, ^  the  experimental  arrangement  and 
initial  progress  on  the  single  source  deposition  of  SiC  was  described 
and  compared  with  earlier  dual-source  work.  The  dual  source  work  had 
yielded  a  non-uniform  composition,  but  in  the  region  of  Si:C  =  1:1,  the 
film  exhibited  columnar  growth  with  'WOO  X  grain  diameter  and  WOO  A 
grain  length  as  111  0-SiC  preferred  orientation  grains  at  ±25°  relative 
rotation.  This  region  of  the  film  had  a  refractive  index  of  2.68  and  a 
sheet  resistivity  of  WQ  cm.  The  initial  single-source  depositions 
were  uniform  in  composition,  appeared  stoichiometric  and  strong  111 
preferred  orientation  was  observed  by  x-ray  analysis  of  films  deposited 
at  700°C  on  (111)  Si  substrates  cleaned  in  situ  by  Argon  ion  bombardment. 
The  present  report  describes  the  continuation  of  the  single-source  work 
and  seven  of  the  eighteen  depositions  to  date  (see  Table  4)  have  been 
selected  to  describe  our  progress. 

The  deposition  apparatus  (see  Fig  13)  and  procedure  are  the 
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same  as  reported  previous ly^  and  three  films  have  been  deposited 
with  cosputtered  tin.  This  involves  rotating  the  heated  substrate 
under  the  ion  gun  for  a  surface  clean  before  rotating  it  through  the 
tin  spray  to  a  position  under  the  SiC  source  which  was  peripherally 
located  in  the  tin  spray  region.  The  SiC  source  was  operated  in  the 
RF  mode,  nearly  always  at  600W,  while  the  Sn  source  was  operated  in 
the  DC  mode  at  low  power  (^20W).  The  system  Argon  pressure  was 
reduced  to  4pm  during  deposition.  At  the  completion  of  deposition, 
the  substrate  was  moved  back  under  the  ion  gun  and  a  slow  cool-down 
period  ensued.  The  latest  depositions  (with  the  tin  added)  showed  the 
best  film  structure  and  properties  to  date. 

(a)  Film  Analysis 

A  number  of  thin  films  have  been  deposited  under  various 
conditions  as  listed  in  Table  4.  The  process  column  describes,  in 
order,  the  SiC  target  power,  substrate  temperature,  source- to- substrate 
distance,  ion  gun  voltage,  ion  etching  chamber  pressure  and  tin. target 
power.  Of  the  work  represented  by  Table  4,  the  depositions  #222-#233 
were  made  during  this  reporting  period.  Detailed  analysis  of  represen¬ 
tative  films  (indicated  with  a  star  in  Table  4)  has  been  started  where 
the  choice  was  made  based  upon  the  presence  of  significant  x-ray  111 
8-SiC  reflections. 
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(b)  Film  Stoichiometry 

Film  compositions  were  measured  by  electron  microprobe  and 
all  samples,  except  those  involving  the  tin  source,  had  profiles  that 
were  flat  to  within  the  reproducibility  of  the  technique.  In  Fig.  14, 
the  data  are  presented  for  the  atomic  percent  Si  found  in  the  represen¬ 
tative  thin  films.  A  standard  of  single  crystal  a-SiC  was  used  to 
determine  the  relative  amounts  of  carbon  and  silicon.  As  seen  in  the 
figure,  films  177,  178,  and  228  are  about  52%  Si.  The  two  samples 
made  by  cosputtering  tin,  231  and  233,  are  quite  different.  These  vary 
from  45%  to  49%  Si  with  the  variation  not  being  fully  understood  at  this 
time  (it  may  relate  to  surface  solubility  and  sticking  coefficient 
changes  because  of  an  expected  tin  thickness  variation) .  The  values 
for  Si  content  are  considered  close  enough  to  50%  to  be  stoichiometric 
SiC,  except  for  deposition  #184  at  200  watts.  In  this  case,  the  Si  is 
^ 45  At%  indicating  a  change  in  the  Si  sticking  coefficient  for  this 
lower  SiC  target  power.  As  shown  later,  some  of  these  samples 
(particularly  178  and  233)  have  significant  B-SiC  111  x-ray  reflections. 

(c)  Refractive  Index 

The  refractive  index  of  the  selected  samples  was  increased 

using  an  automated  Guertner  ellipsoraeter  with  a  HeNe  laser.  To  enhance 

the  variation  in  i|i  for  an  expected  index  ^2.7,  the  measuring  angle  was 

set  to  50°  instead  of  at  the  usual  70°  angle  of  incidence.  These 

results  are  listed  in  Table  5  where  we  see  that  the  refractive  index 

(3) 

values  differ  from  the  literature  range  of  2.65-  2.70  for  SiC. 


iV 


The 


early  depositions,  177  and  178,  are  quite  high  while  the  more  recent 
depositions,  184-233,  are  in  the  proper  range  and  the  low  index  value 
for  sample  180  (on  sapphire)  is  attributed  to  substrate  surface 
roughness. 

(d)  Film  Structure  and  Morphology 

The  thin  film  structure  was  determined  from  x-ray 
diffractometry.  Typically,  the  only  x-ray  reflection  observed  was  the 
8-SiC  111  reflection,  indicating  that  the  films  had  a  strong  texture. 
Accordingly,  Read  camera  photographs  were  taken  in  a  cylindrical  camera 
using  a  15°  angle  of  incidence.  From  conventional  diffractometry,  peak 
heights  were  measured  and  used  to  compare  the  thin  films  to  one  another 
as  shown  in  Table  6.  The  two  films  of  particular  interest  are  178  and 
233,  both  exhibiting  very  strong  0-SiC  111  reflections.  The  angular 
spread  of  the  111  reflection  can  be  measured  on  the  Read  patterns;  i.e., 
the  arc  length  is  approximately  equal  to  the  angle  between  the  111 
reflection  and  the  thin  film  surface  normal.  In  the  Read  photographs, 
other  reflections  were  observed  but  they  were  very  faint  and  Table  6 
presents  a  summary  of  this  analysis. 

The  latest  deposition,  #233,  exhibits  the  highest  degree  of 
preferred  orientation.  For  this  sample,  the  primary  processing  differ¬ 
ences  were  (a)  high  ion-gun  cleaning  voltage  and  (b)  low  power  tin  co¬ 
sputtering.  More  samples  need  to  be  generated  under  these  conditions  to 


confirm  the  importance  of  the  tin  additive. 

For  all  depositions  with  a  substrate  T > 700°C  and  the  substrate 
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surface  had  been  ion- gun  etched,  a  surface  roughening  problem  appeared, 
primarily  along  the  line  of  the  heater.  This  roughening  appears  to  the 
eye  as  a  whitish  area  on  the  film.  We  plan  to  change  heaters,  from  the 
quartz  lamp  type  to  a  carbon-cloth  type  for  the  next  few  depositions  in 
order  to  gain  a  better  perspective  on  this  phenomenon. 

All  of  our  samples  are  being  analyzed  by  TEM  for  orientation  and 
grain  size  and  by  SEM  for  a  detailed  look  at  the  surface.  The  TEM 
selected  area  diffraction  patterns  (SADP's)  will  provide  further 
evaluation  of  the  crystalline  perfection.  To  date,  samples  177  and 
184  have  undergone  TEM  analysis  and  the  SADP  for  the  interfacial  region 
of  177  reveals  a  high  degree  of  crystallinity, 

although  it  has  some  curious  features.  A  two- layer  film  is  suggested 
by  the  overlap  of  the  001  and  112  type  zone  axes,  the  latter  of  which 
is  twinned.  The  dark  field  image  revealed  small  grains  'WOO  X  diameter 

O 

and  larger  grains  WOO  A  diameter.  The  111  type  zone  axis  was  viewed 
for  #184  and  the  spot  pattern  suggests  less  perfection  than  for  #177. 

(e)  Sheet  Resistivity 

Of  the  seven  thin  films  chosen  for  detailed  study,  only  #180 
is  on  sapphire  so  the  sheet  resistivity,  pg,  could  be  measured  with  a 
four-point  probe.  For  the  thin  films  on  silicon,  the  substrate  is  more 
conducting  than  the  thin  film  so  the  electrical  measurements  are  more 
ambiguous.  For  this  section,  all  the  single  source  films  deposited  on 
sapphire  have  been  measured  as  have  the  six  selected  depositions  on  Si. 
There  is  a  substantial  difference  in  pg  depending  on  the  substrate  as 
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expected.  Hie  values  for  SIC  on  Si  are  far  too  low  to  represent  the 
SiC  (typically  'V/ 0.1ft  cm  for  a  lym  film).  This  suggests  that  the  region 
sampled  included  a  significant  depth  into  the  silicon  substrate.  In 
contrast,  for  depositions  on  sapphire,  exhibits  reasonable  values. 

A  listing  of  the  pg  values  for  thin  films  on  sapphire  substrates  is 
given  in  Table  8.  The  expected  order  of  resistivity  for  an  undoped 
film  of  about  lym  thickness  should  be  'v-lO^ft  cm.  ^  Two  of  the  thin 
films,  #122  and  #124,  are  close  to  this  range.  The  lower  pg  value  for 
#125  reflects  the  Si  deficiency  of  this  film  deposited  at  200W. 

(f)  Conclusions  and  Future  Steps 

The  SiC  deposition  process  is  continually  improving. 
Subsequent  depositions  will  place  primary  emphasis  on  the  tin  source 
parameters,  on  improved  surface  cleaning  technique  and  higher  substrate 
temperatures.  In  addition  to  continuing  the  TEM  and  starting  the  SEM 
work  on  the  samples  discussed  in  this  report,  optical  transmission 
measurements  are  underway  for  the  films  on  sapphire.  Optical  .trans¬ 
mission  may  be  possible  for  some  of  the  films  on  Si  by  using  a  special 
sample  preparation  technique  (similar  to  TEM  sample  preparation). 

3.  Synthesis  of  (Pd,  Pt,  Ir)  Silicide  Films 

It  is  our  intention  to  consider  the  important  compounds  in  the 
quaternary  system  (Pd,  Pt,  Ir,  Si)  and  fabricate  high  quality  Schottky 
barriers  of  these  selected  silicides.  Ultimately,  using  these  special 
compounds,  we  propose  to  make  multilayer  structures  with  Si  for  IR 


detector  applications.  The  key  steps  in  the  effort  sequence  are  as 


follows: 


(i)  using  multiple  sources,  determine  the  phase 
spreads  for  all  the  binaries  Pd-Si,  Pt-Si, 

Ir-Si  and  pseudo  binaries  (Pd,Pt)-Si,  (Pd,Ir)-Si 
and  (Pt,Ir)-Si, 

(ii)  development  of  single  source  deposition  for 
special  compounds  in  the  set, 

(iii)  development  of  grain  size  control  in  films  of 
these  special  compounds, 

(iv)  determination  of  Schottky  barrier  heights  for 
these  Si/ (Pd,Pt, Ir)Si  hetero junctions,  and 

(v)  development  of  multilayer  structures  with  the 
heterojunctions . 


(a)  (Pd,  Pt)  Silicide  Film  Synthesis 

A  literature  survey  has  been  conducted  to  gain  structural, 
phase  equilibrium  and  electrical  characteristics  information  concerning 
these  films. 

The  Pt^Si  infrared  detector  is  a  photoconductor  with  a  high 
dark  current.  At  the  hetero junction  between  a  silicide  and  silicon,  the 
conduction  band  of  Si  bends  up  a  distance  <p  where  $(Pt2Si)  *  0.27ev 
and  $(Pd2Si)  *  0.72ev.  The  absorbtion  of  infrared  photons  by  the 
silicide  at  this  junction  causes  changes  in  conductivity  (X  *4.6|im  for 
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Pt2Si  and  Xc  *  1.7ym  for  Pd2Si) .  Devices  are  usually  made  on  20-50Q  cm 
n-type  (100)  Si  although  electrical  measurements  are  easier  on  2-68  cm 
material  and  better  interface  definition  and  uniformity  are  possible  on 
<111>  Si.  It  is  the  extreme  uniformity  of  large  staring  arrays  that 
make  these  silicides  useful  as  IR  detectors.  One  would  like  to  extend 
this  capacity  to  lower  barrier  height  hetero junctions  to  sample  the 
3-6ym  and  the  8-12ym  ranges. 

The  literature  shows  that  Pd2Si  has  been  grown  epitaxially  on 
Si  and  it  may  be  possible  to  grow  Si  epitaxially  on  Pd2Si.  This  is 

O 

because  the  Si-Si  spacing  in  the  Pd2Si  basal  plane  is  3.75  A  which  is 
close  to  the  Si-Si  spacing  of  3.84  £  in  the  (111)  plane  of  Si.  The 
major  compounds  in  the  Pd-Si  system  are  PdSi  (orthorhombic),  Pd2Si 
(hexagonal),  Pd^Si  (orthorhombic-cementite  structure),  Pd4Si2  (orthorhombic- 
metastable)  and  PdgSi  which  melts  congruently.  PdSi  is  not  an  equilibrium 
phase  below  600°C  and  Pd^Si  is  a  common  amorphous  composition.  Pd2Si  is 
not  a  line  compound  and  there  is  some  evidence  that  it  has  a  hexagonal 
superlattice  with  cell  dimensions  2a  and  8c  (similar  to  Fe2P).. 

Although  one  may  seek  a  hetero junction  between  Si  and  the 
different  silicides  to  obtain  a  lower  barrier  height,  other  mechanisms 
of  barrier  height  modification  exist.  Three  examples  of  this  are 
(1)  ion  bombardment  to  produce  deep  levels  which  lead  to  interface 
states  and  thus  a  recombination  current  barrier  lowering;  (2)  Pd2Si 
growth  into  Si  pushes  impurities  leading  to  a  sharply  peaked  distribu¬ 
tion  at  the  interface  -*■  barrier  lowering,  and  (3)  composition  variations 
-*■  barrier  lowering.  Thus,  the  barrier  height  at  the  Pd2Si/Si  junction 
may  be  due  to  interface  nonstoichiometry. 


On  the  experimental  side,  both  the  Pd  and  Si  sources  have  been 
characterized  with  respect  to  deposition  rate  versus  source  power  and 
the  dual  source  Pd-Si  phase  spread  synthesis.  Figure  16  illustrates 
the  dual  source  sputtering  configuration  while  Figs.  17  and  18  provide 
deposition  rate  data  as  a  function  of  electrical  characteristics  for 
Si  and  Pd  respectively.  In  Fig.  19,  the  dependence  of  relative  Si 
deposition  rate  with  Argon  pressure  and  position  is  presented. 

Combining  this  data,  the  composition  spread,  and  thus  the  phase  spread, 
relative  to  distance  from  the  Pd  source  is  given  in  Fig.  20.  The  pre¬ 
dicted  position  of  occurrence  of  the  Pd2Si  composition  is  also  indicated 
in  Fig.  20  (6  cm  from  the  Pd  source). 

4.  Electrical  Diagnostics  of  SiO^  Films 
(a)  I-V  and  C-V  Characteristics 

The  current  vs  voltage  plots  (I-V)  have  been  made  using  a 
DC  sweep  voltage  while  the  capacitance  vs  voltage  plots  (C-V)  have  been 
made  with  a  DC  bias  sweep  voltage  plus  a  lfttiz,  io  mV  signal.  The  pro¬ 
cedures  for  these  plots  are  extensively  developed  and  well  documented 
for  SiC^.  The  two  important  macroscopic  electrical  parameters  obtained 
from  these  plots  are  the  bulk  resistivity,  p,  from  the  I-V  plot  and  the 
bulk  permittivity,  e,  from  the  C-V  plot.  These  standard  tests  have 

been  made  for  the  samples  listed  in  Table  9  at  a  number  of  different 

(p-type  substrates). 

locations  on  the  sampled  Such  measurements  were  made  both  before 
annealing  (BA)  and  after  annealing  (AA)  at  435°C  for  25  minutes.  The 
voltage  range  was  ±20V  and  several  I-V  measurements  were  made  per  dot 
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until  consistent  readings  were  recorded.  This  stable  measurement 
endured  in  subsequent  testing  of  the  same  dot.  Table  9  gives  a  listing 
of  the  resistance  and  capacitance  measurements  in  the  BA  and  AA  condi¬ 
tion  (here,  N  •  R  =  not  reliable  and  S  •  R  =  same  range  as  the  BA  condition) . 
For  this  type  of  diode,  it  is  the  negative  part  of  the  voltage  sweep  that 
leads  to  current  flow  (negative)  in  the  forward  bias  direction  vis-a-vis 
the  p/n  junction  diode  model.  They  may  be  visualized  as  being  inverted 
through  the  origin  for  comparison  purposes.  gl0  j 

In  Figs.  21a  and  21b,  we  see  the  BA  and  AA  for  X  ~  2/on  dot  #2. 
After  annealing,  only  a  small  change  is  noted  in  the  onset  of  the 
forward  current  but  a  very  large  change  is  seen  in  the  reverse  current 
which  is  the  feature  expected  for  short  circuit  conduction  paths  through 
the  transition  region  (see  Fig.  38d- theory).  Figs.  22a,  22b,  and  22c 
exhibit  the  BA  results  for  other  locations  on  the  wafer  and  illustrate 
the  same  behavior  as  Fig.  21a  (note  the  change  in  scale  for  Fig.  22b). 

This  indicates  electrical  uniformity  of  the  film  and  a  junction  that 
produces  Igen  ■  86nA  ±  10%.  Annealing  of  the  wafer  led  to  a  reverse  cur¬ 
rent  slope  for  dot  #1  almost  identical  with  that  of  Fig.  21b.  In  Figs.  23a 
and  23b,  the  C-V  results  on  this  sample  for  the  BA  and  AA  respectively 
are  given.  Hie  BA  result  reveals  what  appears  to  be  the  normal  junction 
capacitance  plus  two  significant  interface  states,  one  appearing  at 
V  'v  0  and  the  other  at  V^-7  volts.  The  AA  result  also  reveals  both  of 
these  features  with  the  state  at  -7V  increasing  greatly  in  magnitude  at 
the  expense  of  the  depletion  layer  capacitance. 


At  X  =  1.5  for  sample  81-083/  a  similar  behavior  is  exhibited 
except  that  there  is  a  greater  variance  from  dot  to  dot.  Figs.  24a  and 
24b  give  the  BA  and  AA  respectively  for  dot  #1  on  this  sample  and  we 
notice  a  significant  change  in  the  contact  potential  effect  for  this 
sample  compared  to  that  represented  by  Fig.  21.  The  point  of  transi¬ 
tion  to  the  forward  current  is  significantly  changed  in  the  AA  case 
compared  to  the  BA  case  and  is  probably  due  to  the  short  circuiting 
effect  of  the  transition  layer.  For  this  sample,  the  variance  in  the 
key  parameters,  (i)  slope  of  reverse  characteristic,  (ii)  slope  of 
forward  characteristic,  (iii)  onset  voltage  for  the  forward  character¬ 
istic,  etc.,  has  increase  to  ±35%.  The  C-V  plots  in  the  BA  and  AA 
condition  for  this  dot  are  presented  in  Figs.  24c  and  24d  and, 
although  these  samples  do  not  clearly  reveal  special  interface  states, 
other  dots  on  this  sample  do  reveal  such  states. 

As  X  decreases  further,  the  onset  voltage  for  the  transition 
to  the  forward  characteristic  increases  further  and  the  variance  of  the 
key  characteristics,  dot  to  dot  on  the  same  sample,  increases.  In  fact, 
for  X  'v-  0. 8,  the  dot  to  dot  variance  is  greater  than  the  average  variance 
between  two  samples  of  different  X. 

(b)  AC  Conductance  Characteristics 

For  reliable  measurements,  a  finely  adjustable  probe  was 
used  to  contact  the  A1  dots  on  the  oxide  side  while  the  Si  ride  had  a 
large  area  A1  contact  held  to  a  vacuum  chuck  and  the  entire  system  was 


contained  in  a  shielded  metal  box.  Fig.  25  is  the  result  for  sample  81-084  (S10, 
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with  data  analysis  leading  to  C  *  132pF  and  p  =  1.68  x  1080  cm  for 

an  estimated  R  of  280K&.  The  measurement  indicated  by  *  was  made 

after  turning  off  the  vacuum  chuck  that  held  the  sample  to  the  contact 

base  and  the  result  shows  the  addition  of  a  large  series  capacitance 

due  to  a  poor  contact  to  the  base.  Fig.  26,  made  with  a  1  volt  signal 

strength  reveals  the  presence  of  two  coupled  circuits.  The  data  can 

be  described  as  the  superposition  of  two  semicircles;  each  with  its 

center  slightly  below  the  abscissa.  This  indicates  the  existence  of  two 

parallel  RC  circuits  connected  in  series  with  each  circuit  having 

frequency  dependent  parameters.  Assuming  frequency  independent 
(a  reasonable  first  order  approximation), 
parameters/ we  find  =  79pF,  ^  =  50k0  and  Px  =  0.33  x  108fi  cm 

plus  C2  =  252pF,  R2  *  300kfi  and  p2  =  2.0  x  1080  cm. 

In  Figs.  27  and  28,  the  p  and  e  values  for  these  films  have 

been  plotted  versus  X.  We  find  a  considerable  difference  between 

samples  measured  at  a  50mV  vs  a  lOOOmV  signal  strength  as  expected 

if  interface  states  are  strongly  involved  in  the  measurement.  We 
likewise  find  a  significant  difference  between  SiOx  samples,  made  by 
cyclic  layer  deposition  (moving  substrate)  vs  simultaneous  Si/C>2 
deposition  (fixed  substrate) . 


B.  THEORETICAL  RESULTS 

1.  Review  of  the  Electrical  Characteristics 
of  MIX  and  Hetero junction  Diodes' 

It  is  our  intention  to  gain  special  insights  into  the  detailed 
structural  features  of  bulk  films,  semiconducting,  insulating  and 
metallic,  and  into  the  electronic  states  at  the  interface  between 


substrate  and  film  by  careful  electrical  characterization  of  the 
composite  system.  The  main  experimental  techniques  to  be  used  are 
I-V  and  C-V  characterization  of  the  diode  characteristics  followed 
by  AC  conductance  characterization.  In  many  of  the  samples  to  be 
studied,  a  thin  native  oxide  exists  on  top  of  the  substrate  prior  to 
film  deposition.  In  many  cases,  this  has  been  removed  by  Argon-ion 
etching. 

A  useful  analogy  with  which  to  begin  discussion  is  the  Si 
p/n  junction  diode.  The  next  step  toward  reality  is  the  hetero junction 
diode  and  the  most  complex  case  is  the  MIS  diode.  We  shall  proceed 
via  these  steps  and  follow  the  discussion  of  Sze.^  in  Fig.  29,  the 
energy  bands,  Fermi  energy,  Ep,  and  electrostatic  macropotential,  V,  is 
given  for  the  isolated  p-type  and  n-type  material  as  well  as  for  the 
junction  where  electron  spill-over  from  n  to  p  has  produced  a  shift  in 
the  electrostatic  potential  so  that  the  Fermi  energies  become  aligned 
and  the  energy  bands  become  shifted. 

In  Fig.  30a,  the  energy  band  diagram  for  two  isolated  pieces 
of  semiconductor  is  given.  The  two  semiconductors  were  assumed  to 
have  different  band  gaps,  E^,  different  permittivities,  e,  different 
work  functions,  $n,  and  different  electron  affinities,  x-  Work  func¬ 
tions  and  electron  affinity  are  defined  as  that  energy  required  to 
remove  an  electron  from  the  Fermi  level,  Ep,  and  from  the  bottom  of 
the  conduction  band,  Ec,  to  a  position  just  outside  the  material 
(vacuum  level).  When  an  ideal  abrupt  hetero junction  without  interface 
traps  is  found  between  these  semiconductors,  the  energy  band  profile 
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at  equilibrium  is  as  shown  in  Fig.  30b  for  an  n-p  anisotype  hetero¬ 
junction.  Since  the  Fermi  level  must  coincide  on  both  sides  in 
equilibrium,  the  discontinuity  in  conduction  band  edges,  AE  ,  and 

v 

valence  band  edges,  AEy,  is  invariant  with  doping  in  those  cases  where 

E  and  x  are  not  functions  of  doping  (i.e.,  nondegenerate  semiconductors). 
© 

The  energy  band  diagram  for  an  ideal  MIS  structure,  for  V  =  0, 
is  shown  in  Fig.  31b  for  a  p-type  semiconductor.  An  ideal  MIS  diode 
is  defined  as  follows:  (1)  at  zero  applied  bias  there  is  no  energy 
difference  between  the  metal  work  function  <|>n  and  the  semiconductor 
work  function  (i.e.,  equal  Fermi  levels)  so  that  a  flat  band  condition 
obtains;  (2)  the  only  charges  that  can  exist  in  the  structure  under 
any  biasing  conditions  are  those  in  the  semiconductor  and  those  with 
the  equal  but  opposite  sign  on  the  metal  surface  adjacent  to  the 
insulator;  and  (3)  there  is  no  carrier  transport  through  the  insulator 
under  DC  biasing  conditions  or  the  resistivity  of  the  insulator  is 
infinity.  When  an  ideal  MIS  diode  is  biased  with  negative  or  positive 
voltages,  the  three  cases  that  may  exist  at  a  p-type  semiconductor 
interface  for  increasing  voltage  are  (a)  accumulation,  (b)  depletion, 
and  (c)  inversion.  In  Fig.  31c  the  band  diagram  for  the  case  of 
inversion  is  illustrated  where  the  surface  potential,  <J>g,  is  positive 
and  larger  than  the  midgap  value,  ^  (potential  difference  between  the 
Fermi  level  and  the  intrinsic  level  of  the  bulk) .  The  typical  variation 
of  the  space-charge  density  per  unit  area  of  interface  required  to 
produce  this  potential,  Qg,  as  a  function  of  the  surface  potential,  tj>s 
is  shown  in  Fig.  32  for  p-type  silicon  with  = 4  x 10^  cm*"*  at  room 


temperature. 
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The  state  of  charge  and  electric  field  in  the  vicinity  of  the 
p/n  junction  and  the  hetero junction  of  Fig.  30  are  illustrated  in 
Fig.  33.  The  chief  difference  here,  which  has  not  been  illustrated, 
is  that  the  hetero junction  may  contain  a  population  of  charged  inter¬ 
face  states  (traps) .  In  case  (a) ,  we  have  that  the  negative  of  the 
contact  potential  is  simply  the  area  under  the  £(x)  vs  x  triangle 
leading  to 
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The  relative  voltage  supported  in  each  semiconductor  is 


NA2e2 

ND1e1 


(7) 


where  V  ■  +  V2.  Itis  apparent  that  Eqs.  5-7  will  reduce  to  the 

expressions  for  the  p/n  junction  (homo junction)  when  both  sides  of  the 
hetero junction  have  the  same  materials. 

Figure  34a  shows  the  band  diagram  of  an  ideal  MIS  structure. 
The  charge  distribution  is  shown  in  Fig.  34b.  For  charge  neutrality 
of  the  system  it  is  required  that 


^  +  *NAW  "  Qs  W 

where  QM  is  the  charges  per  unit  area  on  the  metal,  is  the  electrons 
per  unit  area  in  the  inversion  region,  qNAW  is  the  ionized  acceptors 
per  unit  area  in  the  space  charge  region  with  space  charge  width  W  and 
Qs  is  the  total  charges  per  unit  area  in  the  semiconductor.  The 
electric  field  and  potential,  as  obtained  by  first  and  second  integra¬ 
tions  of  Poisson’s  equation,  are  shown  in  Figs.  34c  and  34d  respectively. 
In  the  absence  of  any  work  function  differences,  the  applied  voltage 
will  appear  partly  across  the  insulator  and  partly  across  the  semi¬ 
conductors  ;  i .  e . , 


v  -  vi  + 
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where  is  the  potential  across  the  insulator  and  is  given  by 
(see  Fig.  34c) 


Vi  s 


Qsd 


c. 

1 


(10) 


where  is  the  insulator  capacitance. 

Figure  35  shows  the  band  diagram,  the  charge  distributions, 
the  electric  field  and  the  potential  for  an  MIS  diode  with  both  surface 
states  and  insulator  space  charge.  Comparing  this  figure  with  Fig.  34 
for  the  ideal  MIS  diode  we  note  that,  for  the  same  surface  potential 
i Ps>  the  applied  voltage  V  is  reduced  indicating  a  voltage  shift  of  the 
C-V  curve  toward  negative  voltage.  The  voltage  shift,  AV,  due  to  the 
space  charge  density  Pi(x)  is  given  by 
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(a)  I-V  Characteristics 

If  we  now  turn  our  attention  to  the  steady- state  conditions  for 
both  forward  and  reversed  biased  junctions  we  notice  that,  for  the  p/n 
junction  case,  the  change  in  the  contact  potential  difference  and  in  the 
energy  band  picture  is  as  given  in  Fig.  36.  The  transition  width,  W, 
decreases  under  forward  bias  and  increases  under  reverse  bias  and  the 
contact  potential  difference  changes  in  the  same  way.  This  gives  rise 
to  the  classical  I-V  characteristic  for  a  p/n  junction  shown  in  Fig.  37 


wherein  the  total  current  is  given  by 

I  =  |l(gen)|(eqV/kT-  1)  .  (12) 

In  Eq.  (12),  the  applied  voltage  V  can  be  positive  or  negative, 

V  =  Vj  or  V  *  -  V  .  The  current  thus  increases  exponentially  with  for¬ 
ward  bias.  When  V  is  negative,  the  exponential  approaches  zero  and  the 
current  is  given  by  - | I (gen) |  which  is  in  the n to  p  (negative)  direction. 
This  current  is  called  the  generation  current  or  the  reverse  saturation 
current  and  is  due  to  the  resulting  drift  of  minority  carriers,  produced 
by  the  thermal  excitation  of  electron-hole  pairs  within  a  diffusion 
length  of  the  junction,  across  the  junction;  i.e., 

D  D 

-  1 1  (gen)  |  *  -qA(^  Pn  ♦  l5-  n  )  (13) 

p  n  v 

where  A  is  the  cross-sectional  area  of  the  junction,  L  and  D  are  the 

n  n 

diffusion  length  and  diffusivity  of  electrons  on  the  n-side  of  the 
junction  while  n  is  the  minority  carrier  population  on  the  p-side  of 
the  interface.  The  reverse  breakdown  current  is  ascribed  to  either 
(a)  Zener  breakdown  due  to  carrier  tunnelling  across  a  heavily  doped 
junction  under  reverse  bias  conditions,  or  (b)  avalance  breakdown  due 
to  impact  ionization  for  lightly  doped  junctions. 

In  the  foregoing,  it  was  implicitly  assumed  that  all  of  the 
applied  bias  voltage  V  appears  across  the  transition  region  rather 
than  in  the  neutral  n  and  p  regions.  Of  course,  there  will  be  some 
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voltage  drop  in  the  neutral  material  if  a  current  flows  through  it; 
however,  for  moderate  to  heavy  doping  the  resistance  of  this  neutral 
material  is  small  and  the  IR  drop  can  be  neglected.  For  the  MIS  diode, 
on  the  other  hand,  this  will  not  be  true  so  that  the  current  through 
the  semiconductor/insulator  junction  must  be  given  by 

I  *  1 1 (gen) |  [exp (q{V-IR}/kT)  -  1]  (14) 

and  the  qualitative  change  in  the  I-V  characteristic  for  this  case  is 
illustrated  in  Fig.  38a.  We  note  a  bending  over  of  the  current  at  a 
large  forward  bias  as  the  characteristic  associated  with  this  effect. 

We  should  thus  expect  this  behavior  for  thick  oxides  and  high 
resistivity  oxides. 

In  analyzing  the  p/n  junction,  it  is  generally  assumed  that 
recombination  and  thermal  generation  of  carriers  occur  primarily  in 
the  neutral  p  and  n  regions,  outside  the  transition  region.  In  this 
model,  forward  current  of  the  diode  is  carried  by  the  recombination 
of  excess  minority  carriers  injected  into  each  neutral  region  by  the 
junction.  Similarly,  the  reverse  saturation  current  is  due  to  the 
thermal  generation  of  electron-hole  pairs  in  the  neutral  regions  and 
the  subsequent  diffusion  of  the  generated  minority  carriers  to  the 
transition  region,  where  they  are  swept  to  the  other  side  by  the  field. 
In  many  cases  this  model  is  adequate;  however,  for  the  MIS  diode 
junction  a  more  complete  description  may  be  necessary  which  includes 
recombination  and  generation  within  the  transition  region  itself. 


This  effect  produces  a  negative  slope  to  the  reverse  saturation 
current  as  illustrated  in  Fig.  38b. 

Another  important  diode  characteristic  is  noted  in  Fig.  38c 
for  heavily  doped  p-n  junction  diodes  having  different  band  gaps.  The 
I-V  characteristics  appear  more  square  than  exponential,  the  current 
being  very  small  until  a  critical  forward  bias  is  reached  and  then  the 
current  increases  rapidly.  This  occurs  because  of  the  difference 
between  low-level  and  high-level  carrier  injection  in  a  junction. 

When  the  injected  excess  carriers  are  negligible  compared  to  the  equi¬ 
librium  majority  carrier  densities,  the  simple  diode  equation  is  valid. 
However,  the  resulting  current  may  be  very  small  particularly  in  heavily 
doped  junctions  and  at  low  temperatures.  As  the  injection  level  is 
increased  to  the  point  where  the  majority  carrier  densities  deviate 
significantly  from  their  equilibrium  values,  a  more  complete  expression 
must  be  used;  i.e., 

.fLiWW),} 

n 

n 

(15) 

As  the  applied  voltage  approaches  the  contact  potential,  the  current 
increases  rapidly.  The  contact  potential  limits  the  forward  bias  which 
can  be  applied  to  the  junction.  In  effect,  the  contact  potential  serves 
as  a  barrier  to  the  net  injection  of  excess  carriers  across  a  junction. 
If  Vq  is  small,  the  diode  current  can  be  significant  at  fairly  low 
voltage;  however,  if  VQ  is  essentially  the  full  band  gap  voltage,  the 


current  is  negligible  until  the  forward  bias  approaches  Eg/q. 

For  an  MIS  diode  of  Si/SiOx/M,  as  x  is  decreased  from  2  in 
the  SiOx  film,  the  film  becomes  more  conducting,  its  Fermi  energy 
increases,  its  contact  potential  difference  w.r.t.  the  p-Si  substrate 
increases  and  it  contributes  electrons  to  the  substrate  leading  to  a 
greatly  increased  Vq  compared  to  that  illustrated  for  the  homojunction. 
As  x  decreases,  we  should  expect  to  see  the  knee  of  the  forward 
characteristics  on  the  I-V  plot  move  to  larger  values  of  V. 

The  effect  of  a  graded  rather  than  an  abrupt  p/n  junction 
does  not  change  the  qualitative  features  of  our  conclusions.  However, 
it  will  alter  the  quantitative  conclusions.  The  effect  of  short 
circuit  conduction  paths  through  the  transition  region  on  the  SiOx 
side,  on  the  other  hand,  will  produce  qualitative  changes.  In  the 
forward  direction,  this  should  act  like  an  additional  current  that 
is  linearly  proportional  to  the  applied  voltage  and  inversely  pro¬ 
portional  to  the  leakage  resistance.  In  the  reverse  direction,  this 
should  act  like  a  generation  current  that  is  linearly  dependent  on  the 
applied  voltage.  Fig.  38d  illustrates  the  qualitative  type  of 
behavior  expected. 

(b)  C-V  Characteristics 

There  are  basically  two  types  of  capacitance  associated 
with  a  junction:  (1)  the  junction  capacitance  due  tc  the  dipole  in 
the  transition  region  and  (2)  the  charge  storage  capacitance,  or 
diffusion  capacitance,  arising  from  the  lagging  behind  of  the  voltage 
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as  the  current  changes,  due  to  charge  storage  effects.  The  junction 
capacitance  is  dominant  under  reverse  bias  conditions  while  the  charge 
storage  capacitance  is  largest  when  the  junction  is  forward  biased. 

For  a  p/n  junction,  the  junction  capacitance,  C.,  is  given  by 


£A  A  ,  2qe  NDNA  % 
j  "  W  *  2  l(VQ- v)  Nd  +  NaJ 


(16) 


where  v  is  the  AC  voltage,  eq.  4  has  been  used  for  W  and  the  quasi 
steady  state  approximation  has  been  used  to  obtain  the  excess  carrier 
density  associated  with  changes  in  v  with  time. 

For  forward  bias,  the  charge  storage  capacitance,  Cs,  is 
calculated  using  a  small  signal  AC  analysis  of  the  injected  carrier 
distribution  and  one  finds  that 
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In  eqs.  17,  Xp  is  the  average  lifetime  of  a  hole  in  the  n-type 
material.  The  factor  Gg  is  an  AC  conductance.  The  charge  storage 
capacitance  can  be  a  serious  limitation  for  p/n  junctions  at  high 


frequency.  In  Fig.  39a,  the  C-V  plot  for  the  p/n  junction  under 
reverse  bias  is  essentially  the  capacitance  associated  with  the 
depletion  layer. 


The  capacitance  for  the  hetero junction  is  given  by 
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For  the  ideal  MIS  diode,  the  total  capacitance,  C,  of  the  system  is 
the  series  combination  of  the  insulator  capacitance,  Ci(=ei/d)  and 
the  semiconductor  space-charge  capacitance,  CD;  i.e.. 


C  = 


CiCD 

Ci*CD 


(19) 


For  a  given  insulator  thickness,  d,  the  value  of  is  constant  and 
corresponds  to  the  maximum  capacitance  of  the  system.  The  semiconductor 
capacitance,  Cp,  depends  upon  voltage  according  to  the  following 
equation 
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where 


3  -  q/kT,  Ld  -  (Ze^qp^S) 


and 


F(8t|>,  5^)  -  [(e~^  ♦  Bifi  -  1)  +~^-(e^  -  8<J>  -  1)]**  >  0  (21) 


po 


po 


At  the  flat-band  condition,  \ps  -  0  and  CQ  becomes 


CD(flat-band)  ■  ^  £s^LD 


(22) 


The  complete  description  of  C  for  the  ideal  MIS  curve  is  given  in 
Fig.  40  where  the  total  capacitance  at  the  flat-band  condition  is 
given,  for  p-type  Si,  by 


«i/id  +  <r> 

s 


(23) 


On  the  left  side  of  Fig.  40,  we  have  an  accumulation  of  holes 
and  therefore  a  high  differential  capacitance  of  the  semiconductor  so 
that  C  Cj.  The  capacitance  goes  through  a  minimum  in  the  depletion 
region  and  increases  again  as  the  inversion  layer  of  electrons  forms 
at  the  surface.  The  increase  of  that  capacitance  is  dependent  upon 
the  ability  of  the  electron  concentration  to  follow  the  applied  AC 
signal.  This  is  only  possible  at  low  frequencies  where  the  recombi¬ 
nation-generation  rates  of  minority  carriers  (here,  electrons)  can 
keep  up  with  the  small  signal  variation  and  lead  to  charge  exchange 
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in  the  inversion  layer  in  step  with  the  measurement  signal.  As  a 
consequence,  MIS  curves  measured  to  high  frequencies  do  not  show  the 
increase  of  capacitance  on  the  right  side  (Fig.  40b). 

For  the  non- ideal  MIS  diode  (Fig.  35),  analysis  of  the  space- 
charge  region  yields  a  shift,  AV,  in  the  silicon  surface  potential  due 
to  the  different  space  charge  densities,  p^(x),  in  the  insulator  (see 
Eq.  11).  The  distance  x  in  this  formula  is  measured  from  the  metal/ 
oxide  interface  so  it  should  be  obvious  that  charges  near  the  silicon 
interface  will  give  rise  to  the  largest  voltage  shift.  When  p^(x)  is 
independent  of  the  externally  applied  voltage,  the  effect  of  the 
localized  charges  on  the  MIS  C-V  curve  is  a  parallel  shift  of  the 
curve  along  the  voltage  axis  by  the  amount  illustrated  in  Fig.  39a. 
However,  when  a  part  of  p^x)  is  dependent  on  the  applied  voltage, 
as  is  the  case  for  interface  states,  one  has  more  difficulty  because 
the  interface  state  capacitance  must  be  extracted  from  the  measured 
capacitance. 

When  a  voltage  is  applied,  the  interface  levels  will  move  up 
and  down  with  the  valence  and  conduction  bands  while  the  Fermi  level 
remains  fixed.  A  change  of  charge  in  the  surface  state  occurs  when 
it  crosses  the  Fermi  level.  This  change  of  charge  will  contribute  to 
the  MIS  capacitance  and  alter  the  ideal  MIS  curve.  The  basic  electrical 
equivalent  circuit  incorporating  the  surface  state  effect  is  shown  in 
Fig,  41a.  Here,  Cg  and  Rg  are  the  capacitance  and  resistance  associated 
with  the  surface  states  and  are  functions  of  the  surface  potential,  i^g. 
The  product  CgRg  is  defined  as  the  interface  state  lifetime,  x,  and 


determines  the  frequency  behavior  of  the  interface  states.  The 
parallel  branch  of  the  equivalent  circuit  of  Fig.  41a  can  be  converted 
into  a  frequency  dependent  capacitance  in  parallel  with  a  frequency 
dependent  conductance  as  shown  in  Fig.  41b  where 


and 


G  Cut 
_£  »  _i _ 

w  1  ♦  w2t2 


Thus,  the  differentiation  of  the  interface  state  capacitance 

from  the  measured  capacitance  as  in  Fig.  41a  does  not  apply  to  the 

equivalent  parallel  conductance  measurement  because  conductance  arises 

solely  from  the  steady-state  loss  due  to  the  capture  and  emission  of 

carriers  by  interface  states  and  is  thus  a  more  direct  measure  of  these 

properties.  Both  capacitance  and  equivalent  parallel  conductance,  as 

functions  of  voltage  and  frequency,  contain  identical  information  about 

interface  states.  However,  as  Nicollian  and  Goetzberger^  have  shown, 

greater  inaccuracies  arise  in  extracting  this  information  from  a 

capacitance  measurement.  This  is  illustrated  in  Fig.  39b  which  shows 

both  the  measured  capacitance  and  the  equivalent  parallel  conduction 

cn  p-type  Si/SiC^  having  an  'v*  2  x  lO^cm"^  and  an  interface  state 
11  -2 

density  in  the  10  cm  -  eV  range  and  measured  at  5-kHz  and  100  kHz> 

The  largest  capacitance  spread  is  14%  while  the  magnitude  of  the 
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conductance  peak  is  over  one  order  of  magnitude  in  this  frequency  range. 
We  note  that  conductance  measurements  yield  more  accurate  and  reliable 
results  particularly  when  the  density  of  states  is  low  as  in  the 
thermally  oxidized  Si  system.  In  a  system  like  Si/SiO  ,  one  expects 
the  density  of  states  to  be  quite  large  and  there  may  be  a  number  of 
important  states  to  consider.  Thus,  one  might  expect  bumps  to  appear 
in  the  C-V  plots  in  certain  voltage  ranges  for  the  Si/SiOx  system. 

Before  closing  this  section,  we  should  also  note  that  the  fixed 
oxide  charge  reported  in  the  Si/SiC^/M  system  causes  a  voltage  shift, 

AV,  via  E.  11  due  to  the  product  xp^^  where  >  0.  The  identical 
effect  would  occur  for  a  suitable  dipole  layer  with  the  positive  charge 
layer  being  closer  to  the  Si  since  x+p+  -  x"p"  >  0  for  this  case. 

Thus,  this  same  effect  could  be  accounted  for  by  a  thin  layer  of 
altered  dielectric  constant  at  the  interface  (decreased) .  We  thus  see 
that  in  Fig.  41  might  perhaps  be  thought  of  as  two  capacitors  in 
series.  If  the  oxide  also  contains  mobile  ionic  charges,  then  one  of 
these  capacitances  at  lease  will  exhibit  a  frequency  dependence  at  low 
frequency. 

(c)  AC  Conductance  Characteristics 

The  most  straightforward  expressions  of  useful  electrical 
quantities  are  as  complex  numbers  in  the  frequency  domain.  A  valuable 
technique  for  the  representation  of  these  quantities  is  to  plot  the 
real  versus  the  imaginary  parts  of  the  quantity  at  various  frequencies. 
The  frequency  dispersion  of  these  values  then  produces  geometric  loci 


of  points  which  are  combinations  of,  or  variations  on,  semicircles 
or  vertical  lines.  These  shapes  allow  for  extrapolation  to  critical 
frequency  ranges  at  which  actual  property  values  for  particular  net¬ 
work  components  can  be  obtained.  As  an  example,  consider  Figs.  42a 
and  42b,  where  two  simple  networks  are  shown  along  with  their  corre¬ 
sponding  plots  of  imaginary  (X)  versus  real  (R)  impedance  components. 

Figures  42c  and  42d  illustrate  the  effect  of  adding  one  more 
parameter  to  the  circuit.  Note  that,  in  Fig.  42c,  the  R^Cj  and  Rj-C2 
relaxation  times  are  sufficiently  different  that  the  semicircle  and 
straight  line  are  distinct.  However,  if  Cj  and  C2  are  similar  (say 
within  two  orders  of  magnitude),  the  shapes  will  overlap.  The  analysis 
of  Fig.  42d,  then,  consists  of  separating  out  contributions  to  the 
curve  shape  which  corresponds  to  a  particular  relaxation  time  or 
R-C  combination. 

The  complexity  of  any  real  system  makes  this  a  difficult 
problem  which  is  compounded  by  the  limitation  of  the  experimentally 
available  frequency  range.  Often  not  realized,  however,  is  the 
considerable  complementary  information  obtained  by  converting  impedance 
values  to  those  of  other  electrical  quantities  such  as  admittance  and 
permittivity  (complex  dielectric  constant) .  Since  admittance  is  the 
reciprocal  of  impedance,  a  series  R-C  combination  now  plots  as  a  semi¬ 
circle  in  the  complex  plane  while  a  parallel  combination  gives  a 
straight  line.  The  permittivity  plot  offers  the  additional  advantage 
of  extrapolating  to  values  of  dielectric  constant  on  the  abscissa.  The 
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simple  relations  between  these  quantities  are 


Z*  =  impedance  -  R  -  jX 

(26) 

Y*  =  admittance  =  G  +  jB 

(27) 

e*  =  permittivity  =  e '  - 

je" 

(28) 

2 

G  =  conductance  =  R/(R 

*x2) 

(29) 

2 

B  =  susceptance  =  X/(R 

/—■v 

CN 

X 

+ 

(30) 

e'  =  dielectric  constant 

=  B/weQ 

(31) 

e"  =  dielectric  loss  =  G/coeo 

(32) 

where  R  is  resistance,  X  is  reactance,  u>  is  the  radial  frequency  and 
eQ  is  the  permittivity  of  free  space.  It  will  also  be  seen  that,  over 
a  given  frequency  range,  one  plot  may  prove  much  more  valuable  than 
another  for  the  purposes  of  extrapolation. 

If  all  the  components  of  a  real  equivalent  circuit  were  discrete 
frequency  independent  parameters,  the  analysis  would  be  relatively 
straightforward.  This,  however,  is  never  the  case  and  considerable 
controversy  has  been  generated  over  apparently  frequency-dependent 
resistances  and  capacitances  which  cause  complex  plane  circles  to 
appear  to  be  centered  below  the  real  axis  and  straight  lines  to  be 
inclined  from  the  vertical. 


f  7) 

Warburg *■  1  first  considered  the  electrical  response  of  the 

interface  between  a  metal  electrode  and  a  liquid  electrolyte.  He 

found  that,  at  a  single  frequency,  X  =  R  and  that,  over  a  limited 

-h  f  81 

frequency  range  X,  Rau  \  Although  later  work1  J  showed  this  to  be 
only  approximately  true,  it  is  still  generally  accepted  as  Warburg's 
law  that 

-a  -a 

R  =  AjW  1  ;  X  =  A2u>  (33) 

where  a^,  ^  0.5  and  A^  ^  A-2-  This  type  of  behavior  can  be  caused 

by  Faradaic  processes  in  which  current  crosses  the  interface  by  means 
of  an  electrochemical  reaction.  This  is  in  contrast  to  a  non- Faradaic 
process  in  which  the  charged  particles  do  not  cross  the  interface  but 
instead  charge  and  discharge  the  electrical  double  layer.  Many  other 
examples  of  Faradaic  processes  exist.  1  In  any  event,  considering 
the  frequency  dependent  aspects  of  the  net  capacitance  in  MIS  diodes, 
we  are  well  to  consider  the  use  of  a  diffusional  admittance  of  the  form 

*  al  a2 

Yz  =  Ai  u  1  +  jA2  w  (34) 

The  foregoing  plus  the  voltage  sensitive  circuit  parameters 
of  Fig.  41b  make  the  data  analysis  a  delicate  issue.  A  system  of 
measurement  is  needed  wherein  rapid  representation  of  the  three  key 
complex  plane  plots  and  rapid  discrimination  of  key  circuit  components 
from  these  plots  can  be  obtained.  One  wishes  to  then  rapidly  subtract 
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these  known  contributions  from  the  data  and  then  replot  the  remainder. 

One  then  iterates  the  previous  step  until  all  the  circuit  components 
are  identified  and  quantitatively  evaluated. 

2.  Film  Formation  Modeling  via  Computer  Simulation 

An  important  adjunct  to  the  present  experimental  program  is 
the  use  of  a  large  high  speed  computer  to  simulate  semiconductor  or 
insulator  film  formation  on  a  specific  substrate.  The  rapidly  increasing 
capacity  of  new  computers  makes  them  an  ideal  candidate  to  bridge  the 
gap  between  atomic  level  considerations  and  macroscopic  continuum 
(thermodynamic)  level  considerations  of  surfaces  and  interfaces.  It 
allows  us  to  explore  the  domain  of  experimental  variable  space  for  the 
film  deposition  process  and  identify  potentially  successful  areas  for 
experimentation.  Through  the  Joint  Stanford/NASA  Ames  Surface  and 
Microstructure  Institute,  we  have  extraordinary  access  to  such  a 
computer. 

The  objective  of  this  study  is  the  analysis  of  surface  and 
interface  properties  at  microscopic  levels  and  is  designed  to  furnish 
all  the  available  theoretical  tools  to  combine  with  our  experimental 
efforts  in  the  best  possible  way.  Two  types  of  simulations  have  been 
initiated  (1)  parametric  investigations  and  (2)  specific  case  study. 

The  parametric  investigations  are  made  for  the  purpose  of  analyzing 
the  general  behavior  of  different  systems  comparatively  and  to  classify 
the  important  factors  which  influence  both  the  structure  and  the  energy- 
related  properties  of  surfaces  and  interfaces  for  these  systems.  In 
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this  type  of  study,  a  molecular  dynamics  technique  based  upon 
Lennard-Jones  type  pair  interactions  is  used  to  investigate  time- 
dependent  as  well  as  equilibrium  properties.  The  specific  case  study 
deals  with  systems  containing  Si  and  0  atoms.  In  this  case  a  more 
involved  potential  energy  function  (PEF)  is  employed  and  the  system 
is  simulated  via  a  Monte-Carlo  procedure.  This  furnishes  the 
equilibrium  properties  of  the  system  at  its  interfaces  and  surfaces 
as  well  as  in  the  bulk. 

(a)  Parametric  Investigation 

The  following  progress  has  been  made  during  this  first 
six-month  period: 

(1)  A  molecular  dynamics  program  (MOLDY)  based  on  the 
Nordsieck-gear  algorthim  has  been  prepared.  The 
program  is  able  to  deal  with  the  bulk  as  well  as 
with  the  surface  and  interface  regions  for  systems 
containing  up  to  three  different  atomic  species. 

(2)  A  utility  package  has  been  prepared  for  the  evaluation 
of  the  MOLDY  results.  This  package  consists  of  several 
programs  and  uses  the  data  generated  by  MOLDY  as  an 
input.  It  has  the  following  capabilities: 

(i)  The  radial  distribution  function,  virial,  density 
and  energy  profiles  across  the  interface  can  be 


obtained. 
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(ii)  It  can  perform  time  series  analysis:  correlation 
functions  (density  or  velocity  auto-correlation) 
and  vibrational  frequency  spectra  of  the  particles 
in  motion  can  be  calculated.  The  program  has  the 
flexibility  of  calculating  components  of  a 
vibrational  frequency  spectrum  for  any  number  of 
particles  and  for  any  desired  subsection  of  the 
system.  It  uses  a  fast  fourier  transformation 
routine  and  accepts  velocity  components  of  the 
particles  generated  by  MOLDY.  The  program 
requires  relatively  small  amount  of  CPU  time. 

(3)  A  graphic  software  package  has  also  been  developed  which 
provides  visual  inspection  of  the  whole  system  as  well 
as  of  any  local  region  of  the  system.  It  can  be  employed 
for  generating  a  single  frame  for  any  desired  time  step 
or  it  can  be  used  to  produce  a  movie  for  the  real  time 
simulation  of  individual  atomic  motions. 

The  package  also  includes  a  program  which  can  be 
used  to  generate  trajectories  of  individual  particles 
for  quick  inspection. 

At  the  moment  a  test  case  is  under  investigation  for  the  analysis 
The  system  consists  of  two  different  atomic  species  to  represent  an  :nter 
face  between  a  substrate  and  an  adlayer.  Figure  43  demonstrates  the 
system  schematically.  It  contains  927  particles  (with  381  particles  (Cu) 


in  the  adlayer  and  the  rest,  546  particles  (Ag),  in  the  substrate). 
The  reduced  parameters  used  in  the  molecular  dynamics  procedure 
(representative  of  a  Cu/Ag  system)  were: 


id  -  1-1872 

0*d  =  0.8846 

s*ub  =  1*° 

a*  =  i 

sub 

Time  step,  dt*  =  0.01 

Temperature,  T*  =  0.3 

No  periodic  boundary  condition  has  been  imposed  because  of  the  binary 
nature  of  the  system.  The  particles  (Ag)  within  the  three  bottom  layers 
of  the  substrate  were  held  fixed  during  the  run;  but,  they  were  permitted 
to  interact  with  the  rest  of  the  system  contributing  to  the  total  poten¬ 
tial  energy.  The  test  system  was  run  up  to  1000  time  steps.  The  first 
200  steps  were  used  for  the  temperature  rescaling  process,  and  the  rest 
employed  in  the  analysis.  All  the  thermodynamic  quantities  were  calcu¬ 
lated  considering  only  the  "core"  region  of  the  hexagonal  slab  (to 
eliminate  the  surface  effect).  The  core  included  all  the  atoms  within 
a  radius  of  2.8  a*ub  from  the  center.  Figure  44  demonstrates  mobile 
atoms  during  the  equilibration;  and  Figure  45  is  a  close-up  view  of  the 
same  region  around  the  interface.  Figure  46  shows  trajectories  of  the 
atoms  within  the  top  layer  of  the  substrate,  and  Fig.  47  is  a  trajectory 
plot  for  the  bottom  layer  of  the  adsorbate  (which  is  facing  the  layer 
in  Fig.  46) . 
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Figure  48  demonstrates  the  potential  energy  profile  across 
the  interface  which  is  located  between  the  third  and  the  fourth  peaks. 
The  difference  between  these  peaks  indicates  the  stronger  binding  among 
the  adatoms.  Figure  49  is  the  virial  profile  similar  to  Fig.  48. 
Because  of  its  differential  character,  the  virial  reflects  the  mismatch 
of  the  sites  at  the  interface.  The  difference  in  atomic  sizes  causes 
the  adlayer  side  to  dilate  while  the  substrate  side  contracts.  The 
dilation  is  exhibited  by  a  negative  virial  while  the  contraction  by  an 
added  positive  contribution  as  expected. 


(b)  Simulation  of  Systems  with  Si  and  0  Atoms 

For  a  more  quantitative  representation  of  systems  involving 
Si  and  0  atoms  we  considered  a  general  type  of  potential  energy  function 
which  is  based  on  two-body  and  three-body  interactions.  The  most 
serious  difficulty  involved  in  simulation  studies  of  such  systems  is, 
in  general,  due  to  the  potential  energy  function  and  its  ability  to 
represent  the  system  adequately.  The  analytical  form  of  the  potential 
energy  function  and  the  evaluation  of  its  parameters  that  are  being 
used  are  described  briefly. 

For  a  system  of  N  particles  the  total  potential  energy  i  may  be 
expressed  as: 


where,  ufr^,  f..),  uf?^,  T.,  f^),  ....  u(rit  f ....  T)  are  the  two-, 
three-  dnd  n-body  potentials  respectively; 'r^  denotes  the  position  of 
the  j**1  particle.  Clearly,  the  most  important  term  in  this  expansion 
is  the  first  term  involving  the  two-body  potential.  Therefore,  in  the 
majority  of  the  atomistic  calculations  to  date  only  pairwise  additive 
potentials  have  been  used.  In  this  investigation  not  only  two-body 
but  also  three-body  interactions  are  taken  into  consideration.  Because 
this  expression  (Eq.  35)  will  be  used  in  lengthy  machine  computations, 
the  functional  form  of  u(r“. ,  "?. )  and  u(?. ,  ,  ”£)  must  be  very  simple. 

1  J  1  j  K 

The  two-body  part,  therefore,  is  represented  by  a  Mie-type  potential. 


“<V  *  ss 

1 J  1 J 


with  r^  *  |?i  -"fjl,  while  r°  denotes  the  equilibrium  separation  and  e 
is  the  energy  at  r^  =  r°.  The  exponents  m  and  n  account  for  the 
repulsive  and  attractive  terms,  respectively.  For  the  intermediate  to 
large  separations  the  three-body  interactions  may  be  expressed  as: 


ufr.,  f .,  -*k)  -  E  Zt  •  Gjr.,  *f. ,  -*k)  (37) 

where,  the  summation  includes  all  triple  multiple  interactions  resulting 
from  the  expansion  of  the  third-order  interaction  energy  for  three  atoms. 
Each  term  in  the  summation  is  expressed  as  the  product  of  a  geometrical 
factor  G(r. ,  “r.,  ■?.)  which  depends  on  the  relative  positions  of  the  three 

1  J  K 

atomic  nuclei  and  an  interaction  constant  which  depends  only  on  the 
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atomic  species  involved  in  the  interaction.  Here,  we  employ  only  the 
triple-dipole  interaction  which  has  been  shown  to  be  the  dominant 
contribution: 


(38) 


with 


1 + 3Cos  9. Cos  9. Cos  9, 
_ L  3  k 

trin  -  rik  '  rjk>3 


(39) 


where,  9^,  9^,  9^  and  r„,  r^»  r^  represent  the  angles  and  the  sides 
of  the  triangle  formed  by  the  three  particles  i,  j  and  k. 

Now,  combining  equations  35  through  39,  and  neglecting  the 
four-body  and  higher  terms  one  obtains: 


r2  T3 


(40a) 


with 


,  N  N  c  o  m  o 

1  J  13  13 

i^3 


(40b) 


and 


.  NNN  (1  +  3Cos9.  Cos9 .  Cos9.) 

A  r  r  r  t  ^  J  ^ 


h- Zzz  l  zi 

i  3  k 

Wj/k 


(ri3 


rik 


V 


(40c) 
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In  this  general  form  the  total  potential  energy  $  is  a  function 
of  the  •parameters,  e,  r°,  m,  n  and  Z^,  and  the  atomic  configuration  of 
the  system.  The  parameters  are  assumed  to  be  independent  of  the  atomic 
positions  and  the  geometry  of  the  system.  They  depend  only  on  the  atomic 
properties  of  species  involved. 

For  a  binary  system  of  types  1  and  2  [1  =  silicon,  and 
2  =  oxygen]  with  the  corresponding  number  of  particles  and  N2> 
equation  6  may  be  rewritten  as : 


45 


where. 


2  C^-)” 

j  aj 


r  .  ^  0 
aj 


2  c/-)n 

j  aj 


raj  *  0 


NB  \  r  t  0 

E  l  G(ruj ■  rc.k-  rjk>  {  .  „ 

j  <k  aj 


r  0  represents  the  equilibrium  separation  between  the  two  species  a  and 
6  for  the  two-body  part  (r°g  =  rj^);  "d"  is  a  critical  distance 
parameter  which  is  used  as  a  normalizing  factor.  $  is  independent  of 
the  value  of  4.  roj  -  |?B  -  *j|,  where  ?  denotes  th#  of  „ 

atom  of  type  a  taken  as  the  central  particle  for  the  summation.  The 
total  number  of  particles  in  the  system  is  given  by  N^.Qt  »  Nj  +  N2. 

The  parameters  (such  as  eag,  r°g,  aag,  Zagy  and  TagY)  were 
obtained  using  isolated  microclusters  and  the  bulk  property  data  for 
pure  Si  and  0,  as  well  as  data  for  Si02  forms.  Furthermore,  for  the 
evaluation  of  the  parameters,  the  stability  criterion  was  considered 
that  is  given  as: 


46  - 


Next,  we  used  thermodynamic  data  (namely  cohesive  energies) 
for  a-quartz  and  a-cristabolite.  Employing  equations  6  and  14,  the 
rest  of  the  potential  energy  parameters  were  calculated  as: 

e12  =  42531.7  (°K) 

r°2  =  1.6259  (A) 

Z112  =  5.0086  x  106  (°K  *  A9) 

Z122  =  3.9999  x  106  (°K  *  A9) 

This  set  of  parameters,  furthermore,  used  in  equation  40  and 
the  cohesive  energies  for  the  6-cristabolite  and  6-quartz  were 
calculated.  Results  are  found  to  be  in  good  agreement  with  the 
experiments.  Calculated  energies  were: 

Vcristabalite)  '  '414-9  (ltcll','"le) 

♦(8-quartz)  *  -406'5  <kcal/*ole) 
experimental  cohesive  energies  were  tabulated  in  Table  10. 
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For  further  investigation  of  SiO^  systems,  the  following 
programs  are  under  preparation. 

(1)  Monte-Carlo  Program:  This  program  is  written  to 
handle  systems  containing  up  to  three  different  atomic 
species.  It  uses  equation  40  and  generates  appropriate 
Markov  chain  configurations  of  the  system  in  a  very 
efficient  way.  It  stores  the  list  of  neighbors  using 
the  "large  core  memory"  (in  a  CDC  7600)  that  reduces 
the  computation  time  considerably,  in  particular,  for 
larger  systems.  This  program  is  now  in  its  debugging 
stage. 

(2)  A  Software  Package  for  Structural  Characteristics 
of  SiO^  Systems: 

This  package  contains  basically  two  different 
programs.  The  first  is  to  calculate  the  radial 
distribution  function  of  the  system.  It  provides  the 
radial  distribution  function  associated  with  any  portion 
of  the  system;  and  it  is  able  to  calculate  partial 

(i*e*’  g  (Si-Si) ’  g(Si-0)  °r  g[oV  °r  the  t0tal  g(R) 
of  any  Si02  system.  The  second  program  is  to  calculate 

angular  distributions  for  the  whole  system  or  for  any 

subsection  of  the  system.  In  particular,  it  is  designed 

to  calculate  the  distributions  for  angles  such  as 

((Si  -  0  -  Si)  and  (0  -  Si  -  0)  between  the  neighboring 


atoms. 


49  - 


References 

(1)  W.  A.  Tiller,  "Fabrication  and  Properties  of  Multilayer  Structures," 
Fourth  Semiannual  Technical  Report  for  Period  March  1,  1981- 
August  31,  1981,  ARPA. 

(2)  W.  A.  Tiller  and  T.  W.  Barbee,  Jr.,  "Fabrication  and  Properties  of 
Multilayer  Structures,"  Third  Semiannual  Technical  Report  for 
Period  September  1,  1980-February  28,  1981,  ARPA. 

(3)  Handbook  of  Physics  and  Chemistry,  Chemical  Rubber  Co.,  Cleveland, 
Ohio  197 

(41  Silicon  Carbide,  Eds.  R.  C.  Marshall,  J.  W.  Faust,  Jr.,  and 
C.  E.  Ryan  (U.S.  Carolina  Press,  Columbia,  N.Y.)  1974. 

(5)  S.  M.  Sze,  Physics  of  Semiconductor  Devices  (John  Wiley  §  Sons, 

N.Y. ) ,  1981. 

(6)  E.  H.  Nicollian  and  A.  Goetzberger,  Bell  Syst .  Tech.  J.  46. 

1055  (1967). 

(7)  E.  Warburg,  Annalen  der  Physik  und  Chemie  67.  493  (1899). 

(8)  H.  Fricke,  Physics,  JL  106  (1931). 

(9)  D.  C.  Grahame,  J.  Electrochem  Soc.  99,  370  (1952). 


1 


-  50  - 


Table  1 

Ion  Gun  Etch  Parameters 


1 

Ion  Gun 
to  Substrate 
Distance 
(inches) 

Operating 

Voltage 

(kV) 

Chamber 
Pressure 
(ym  Ar) 

Etching 

Time 

(min.) 

Thickness 

Change 

A 

3 

2.0 

4 

9 

120 

15.6 

' 

3 

2.0 

4 

30 

650 

21.7 

3 

2.0 

8 

15 

320 

18.9 

3 

2.0 

8 

30 

575 

18.7 

4 

0.75 

4 

30 

none 

- 

4 

0.75 

8 

40 

none 

- 

1 

4 

2.0 

4 

20 

245 

12.25 

r 

4 

2.0 

4 

40 

650 

16.2 

4 

2.0 

8 

20 

270 

13.5 

4 

2.0 

8 

40 

580 

14.5 

5 

0.75 

4 

20 

none 

- 

5 

0.75 

4 

30 

none 

- 

5 

0.75 

4 

40 

none 

- 

5 

0.75 

8 

20 

none 

- 

5 

0.75 

8 

40 

none 

- 

— 

5 

2.0 

7.5 

15 

none 

- 

5 

2.0 

7.5 

20 

none 

- 

5 

2.0 

7.25 

40 

130 

3.25 

% 

Note:  (1)  All 

substrates  were  'Wym  thermal  oxide  on  Silicon. 

(2)  Operating  current  was  2  mA 

in  all  cases. 

* 
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Table  2 


Temperature  Dependence  of  Oxide  Depositions 


substrates:  (lll)Si 

source:  3"  Si 
distance:  4" 

ion  etch:  750V,  2  mA,  40  min. 
deposition  power:  500  W  (rf) 
time:  20  min. 

oxygen  pressure:  0.30  ±  0.015  ym 
Argon  pressure:  4.0  ym 


Run  Number 

Temperature 

(°C) 

Thickness 

xb 

Refractive 

Index 

0/Si 
(At  %) 

81-215 

25 

7000 

1.47 

65.6/34.4 

81-219 

25 

7300 

1.46 

66.6/33.4 

81-216 

240 

6550 

1.50 

63.4/36.6 

81-217 

340 

4650 

1.46 

60.1/39.9 

81-218 

500 

3500 

1.53 

55.4/44.6 

81-220 

660 

(2500) 

(1.45) 

57.6/42.4 

81-221 

660 

4500 

1.45 

47.0/53.0 

r 

r 


Table  3 


-  ■ 

Cyclic 

Formation 

of  SiO„  Layers 

s 

Sample 

Number 

Si 

Deposition 
0Rate 
(A  /  sec) 

Rotation 
Period 
(  sec  ) 

A 

Si 

Layer  Number 

Thickness  of 

(A)  Layers 

Exposure 
(torr  sec) 

Gross 

Composition 
Si/0  (At  %) 

i 

81-083 

2.05 

0.52 

1.1 

3461 

1.1  x  10"4 

45:55 

81-084 

2.05 

5.2 

10.7 

156 

1.5  x  10"3 

49.5:50.5 

t- 

81-085 

2.05 

10.2 

20.9 

176.5 

2.5  x  10'3 

57:43 

S 
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Table  4 


Single  Source  Depositions 


Film 

Number 

Substrate 

Process:  SiC  target  power,  T,  source-substrate 
distance,  ion  gun  voltage,  P,  Sn 
target  power. 

122 

A12°3 

600W,  25°C,  3" 

124 

A12°3 

600W,  700°C,  3" 

12S 

A12°3 

200W,  820°C,  3" 

*  177 

(100)Si 

600W,  720°C,  3",  750V,  8y 

*  178 

(lll)Si 

600W,  700°C,  3",  750V,  8y 

*  180 

A12°3 

600W,  660°C,  4",  750V,  8y 

*  184 

(lll)Si 

200W,  700°C,  4”,  750V,  8 y 

210 

(100)Si 

600W,  25°C,  4",  750V,  8y 

211 

(100)Si 

600W,  615°C,  4",  750V,  8y 

212 

(100)Si 

600W,  715°C,  4'*,  750V,  8y 

222 

(lll)Si 

600W,  25°C,  5",  750V,  8y 

223 

(lll)Si 

600W,  750°C,  5",  750V,  8y 

224 

(lll)Si 

600W,  700°C,  5",  750V,  8y 

225 

(lll)Si 

600W,  615°C,  5",  750V,  8y 

*  228 

(lll)Si 

600W,  7508C,  5",  2000V,  4y 

229 

(lll)Si 

600W,  750°C,  5",  1000V,  4y,  145W 

*  231 

(111) Si 

600W,  750°C,  5",  1000V,  4y,  18W 

*  233 

(lll)Si 

600W,  750°C,  5",  4000V,  4y,  16W 
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Table  S 

Refractive  Index  Measurements  of  Selected  SiC  Thin  Films 


Sample 

Number 

Substrate 

Refractive 

Index 

Stoichiometry 
Si  :  C 

177 

(100)Si 

3.03 

52:48 

178 

(lll)Si 

3.21 

52:48 

180 

ribbon 

A12°3 

2.34 

51:49 

184 

(lll)Si 

2.71 

45:55 

228 

Clll D  Si 

2.76 

53:47 

231 

(lll)Si 

2.54 

47:53 

233 

(lll)Si 

2.73 

48:52 

Table  6 


Film 

Number 

177 

178 

180 

184 

228 

231 


Representative  SiC  Depositions  and  Relative  X-ray  Peak  Heights 


Substrate  Process  Relative  X-ray  Peak 


233 


(lll)Si  600W,  750°C,  5”,  4kV,  4y,  18W' 


100 


Table  7 


Read  Camera  Analysis 


arc  length  of  111 

Sample  6-SiC  reflections  present  deviation  from 

Number  111  220  311  200  surface  normal 
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Table  8 

Sheet  Resistivity  of  SiC  Thin  Films  Measured 
_  with  a  4- Point  Probe 


Deposition  SiC  Film  Sheet 

Temperature  Target  Thickness  Resistivity  Si:C  Film 
Sample  Substrate  °C  Power  pm  ft  cm  Composition 


122 

1T02  ai2o3 

25 

600W 

0.85 

265 

51:49 

124 

iTo2  ai2o3 

700 

600W 

0.52 

320-590 

50:50 

125 

1102  A1203 

820 

200W 

0.55 

6-24 

46:54 

180 

ribbon 

A1,0, 

660 

600W 

1.40 

260 

51:49 

2  3 


Table  10 


a-Crist. 

a- quartz 

8-Crist. 

8- quartz 

R(Si-Si) 

3.0767 

3.0579 

3.0880 

3.0982 

R(Si-O) 

1.60 

1.60 

1.55 

1.63 

R(O-O) 

m 

2.60 

2.63 

2.53 

2.60 

all 

4.0575 

4.1148 

4.0388 

4.1048 

m 

a12 

9722.33 

8328.37 

9796.98 

9870.78 

am 

22 

41.3494 

36.7629 

42.0266 

48.1561 

m 

a21 

4861.17 

4164.88 

4898.46 

4929.74 

all 

5.2401 

5.6759 

5.0999 

5.6325 

a?2 

202.442 

188.962 

202.892 

204.702 

a» 

19.5795 

20.2410 

19.3041 

21.5124 

a21 

101.221 

94.4946 

101.444 

102.302 

rill 

5.7840 

8.3242 

4.94088 

8.0688 

T112 

-191.313 

-117.663 

-202.276 

-175.685 

ri22 

174.512 

195.381 

167.158 

195.231 

r222 

62.0504 

76.5891 

57.5318 

73.9876 

r22l 

174.555 

195.427 

167.131 

195.237 

r211 

-47.8405 

-29.4278 

-50.5788 

-43.9041 

Cohesive 

Energy  -412.14  -412.540  -417.43  -417.85 

(kcal/mole) 


Figure  Captions 


1.  (a)  Arrangement  of  deposition  source,  Ar-c leaner,  oxygen  supply  ring 

and  substrate  heater  used  in  the  reactive  sputtering  of  SiOx 
films. 

(b)  Expanded  view  of  the  quarts  ring  for  supplying  02  to  the  rec¬ 
tangular  substrate  surface. 

2.  Ion-gun  etching  profile  for  the  following  conditions: 


Run 

Power 

Distance 

Chamber  ; 

Pressure 

Etch 

Time 

Etch 

Rate 

— 

2kV,  2mA 

3  in 

8  um(ar) 

30  min 

18.7  1/mln 

1 

2kV,  2mA 

i 

3  " 

8  "  " 

!  15  " 

18.9  " 

— —  •  • 

2kV,  2mA 

3  " 

4  »»  <» 

30  ” 

21.7  ” 

2kV,  2mA 

3„ 

;  4  "  ” 

! _ i 

9  "  ! 

1  ! 

1 _ 1 

15.6  " 

3.  Index  of  refraction  profile  for  the  same  conditions  as  in  Fig.  2. 

4.  Configuration  of  reactive  sputtering  deposition  mask. 

5.  SiO^  thickness  vs.  deposition  temperature. 

6.  Topographic  map  of  deposit  thickness  via  elllpsometry  for  samples 
82-030  and  82-031. 

7.  Infrared  transmission  spectra  for  typical  amorphous  S102. 

8.  IR  Transmission  spectra  for  amorphous  Si02  at  several  deposition  temper¬ 
atures  Illustrating  the  wavenumber  shift. 

9.  Thickness  profiles,  determined  by  the  Tally-step  method,  for  SiOx  formed 
via  the  cyclic  deposition/reaction  technique. 

10.  Index  of  refraction  profiles,  determined  ellipsometrically,  for  S10x 
formed  via  the  cyclic  deposition/reaction  technique. 
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11.  Plot  of  ZnO  stoichiometry  profile  over  the  substrate  for  the  simul- 
taaeous  deposition/reaction  technique. 

12.  Plot  of  index  of  refraction  variation  for  the  sample  in  Fig.  11. 

13.  Single-source  SiC  sputtering  system  with  ion  gun  for  substrate  surface 
cleaning  and  with  portable  tin  source  in  a  likely  position. 

14.  Stoichiometry  of  selected  thin  films  as  measured  by  electron  microprobe. 

15.  Figure  withdrawn  from  report. 

16.  Dual  source  sputtering  configuration  used  in  these  studies. 

17.  Silicon  deposition  rate  as  a  function  of  R.F.  power. 

18.  Palladium  deposition  rate  as  a  function  of  current  of  the  D.C.  power 
source. 

19.  Film  deposition  characteristics  for  Si  for  various  argon  pressures. 

20.  Relative  Pd/Si  deposition  rate  as  a  function  of  distance  from  the  Pd 
source. 

21.  1-V  plots  for  SiOx  samples  both  before  annealing  (BA)  and  after  annealing 
(AA). 

22.  I-V  plots  (BA)  for  SiOx  at  different  locations  on  the  same  sample. 

23.  C-V  plots  for  SiOx  both  BA  and  AA  at  one  location. 

24.  I-V  and  C-V  plots  for  SiOx  at  one  location  (BA  and  AA). 

25.  A-C  conductance  plots  (reactance  vs.  resistance)  for  SiO  (BA)  at  an 

A 

oscillator  amplitude  of  50  mV. 

26.  A-C  conductance  plot  for  SiOx  at  an  oscillator  amplitude  of  IV. 

27.  SiOx  resistivity  as  a  function  of  stoichiometry  (BA)  for  both  the  fixed 
substrate  (simultaneous  deposition)  and  the  moving  substrate  (sequential 
deposition)  cases. 
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28.  SiO^  deliectric  constant  (BA)  for  the  fixed  substrate  (F.S.)  and  moving 
substrate  (M.S.)  cases. 

29.  Properties  of  an  equilibrium  p-n  junction:  (a)  isolated,  neutral  regions 
of  p-type  and  n-type  material  and  energy  bands  for  the  isolated  regions; 
(b)  junction,  showing  space  charge  in  the  transition  region  H,  the  re¬ 
sulting  electric  field  £  and  contact  potential  Vq,  and  the  separation  of 
the  energy  bands;  (c)  directions  of  the  four  components  of  particle  flow 
within  the  transition  regions  and  the  resulting  current  directions. 

30.  (a)  Energy-band  diagram  for  two  Isolated  semiconductors  in  which  space- 

charge  neutrality  is  assumed  to  exist  in  each  region. 

(b)  Energy-band  diagram  of  an  ideal  n-p  anisotype  heteroj unction  at 
thermal  equilibrium. 

31.  (a)  Metal- insulator-semiconductor  (MIS)  diode. 

(b)  Energy-based  diagram  of  ideal  MIS  diode  at  V  «  0  for  a  p-type  semi¬ 
conductor. 

(c)  Energy-band  diagram  of  ideal  MIS  diode  at  V  >  0. 

32.  Variation  of  space-charge  density  in  the  semiconductor  as  a  function  of 

15  -3 

the  surface  potential  ip  for  p-type  silicon  with  N  ■  4  x  10  cm  at 

S  A 

room  temperature;  <J;  is  the  potential  difference  between  the  Fermi  level 

D 

and  the  intrinsic  level  of  the  bulk  semiconductor. 

33.  (a)  Space  charge  and  electric  field  distribution  within  the  transition 

region  of  a  p/n  junction  with  N ,  >  N  :  top,  the  transition  region 

d  ft 

with  x  ■  0  defined  at  the  metallurgical  junction;  middle .  charge 
density  within  the  transitional  region,  neglecting  the  free  carriers; 
bottom,  the  electric  field  distribution,  where  the  reference  direc¬ 


tion  for  £  is  arbitrarily  taken  as  the  +x  direction. 
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(b)  The  same  for  Si/SiOx  he teroj unction. 

34.  (a)  Band  diagram  of  an  ideal  MIS  diode 

(b)  Charge  distribution  under  inversion  condition 

(c)  Electric  field  distribution 

(d)  Potential  distribution. 

35.  An  MIS  diode  with  fixed  oxide  charge  and  oxide  trapped  charge. 

(a)  Band  diagram 

(b)  Charge  distribution 

(c)  Field 

(d)  Potential. 

36.  Effects  of  bias  at  a  p-n  junction;  transition  region  width  and  electric 
field,  electrostatic  potential,  energy  band  diagram,  and  particle  flow 
and  current  directions  within  W  for  (a)  equilibrium,  (b)  forward  bias 
and  (c)  reverse  bias. 

37.  I-V  characteristic  of  a  p/n  junction. 

38.  Various  non- ideal  p/n  junction  diode  effects. 

39.  (a)  Ideal  p/n  junction  C-V  plot  plus  real  Si/Sit^  heteroj unction  C-V 

plot  illustrating  the  flat  band  shift,  AV,  due  to  the  charge  effects 
at  the  non- ideal  junction. 

(b)  Capacitance  and  equivalent  parallel  conductance  measurement  at  5  kHz 

and  100  kHz  on  a  p-type  sample  having  an  acceptor  density  of  2.08 

16  -3  11  —2  ”1 
x  10  cm  and  interface  state  density  in  the  10  cm  -eV  range. 

40.  MIS  capacitance-voltage  curves:  (a)  low  frequency,  (b)  high  frequency, 

(c)  deep  depletion. 

41.  (a)  Electrical  equivalent  circuit  of  the  MIS  incorporating  the  surface 

state  effect  (Cc  and  R  are  capacitance  and  resistance  associated 

w  J 
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with  the  surface  states). 

(b)  Conversion  of  parallel  branch  of  (a)  into  a  frequency  dependent 

capacitance  in  parallel  with  a  frequency  dependent  conductance  Gp. 

42.  Complex  impedance  plots  for  (a)  series  R-C  circuit,  (b)  parallel  R-C 
circuit,  (c)  C2  «  Cj,  (d)  C2  ~  C^ 

43.  Schematic  representation  of  atomic  system  for  interaction  between  copper 
and  silver. 

44.  Simulation  of  Cu  movement  on  Ag(lll). 

45.  Close-up  view  of  Fig.  44. 

46.  Trajectories  of  the  atoms  within  the  top  layer  of  the  Ag(lll)  substrate 
at  the  Cu/Ag  interface. 

47.  Trajectories  of  the  atoms  within  the  bottom  layer  of  the  Cu  adsorbate 
at  the  Cu/Ag  interface. 

48.  Potential  energy  profile  across  the  interface  illustrating  the  stronger 
binding  among  the  adatoms. 

49.  Virial  profile  across  the  interface  reflecting  the  mismatch  of  the  sites 
at  the  interface  (adlayer  side  dilates  while  the  substrate  side  con¬ 


tracts)  . 
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FIGURE  40 
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